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A largely model-independent measurement of the inclusive electron momentum spectrum and 
branching fraction for semileptonic decays of B mesons is presented based on data recorded at the 
T(AS) resonance with the BABAR detector. Backgrounds from secondary charm decays are separated 
from prompt B decays using charge and angular correlations between the electron from one B meson 
and a high momentum electron tag from the second B meson. The resulting branching fraction is 
B(B — > Xev) = (10.87 ± 0.18(stat) ± 0.30(syst))%. Based on this measurement we determine the 
CKM matrix element | V^b | - 



PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 
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Measurements of semileptonic B meson decays are a 
good way to determine the CKM matrix elements \V c b\ 
and \V u b\i two of the parameters of the Standard Model. 
For \V c b\, analyses of exclusive and inclusive decays have 
resulted in comparable precision. While most measured 
values of B(B —> Xev) are below 11% [[[], theoretical cal- 
culations including perturbative QCD contributions pre- 
dict values of 12% or above B. 

The measurement presented here employs the method 
introduced by ARGUS @ and later used by CLEO @, 
in which BB events are tagged by the presence of a high 
momentum lepton. As a tag, we choose electrons with 
momentum p* in the interval 1.4 to 2.3 GeV/c, where p* 
is measured in the center-of-mass frame. A second elec- 
tron in the event is taken as the signal lepton for which 
we require p* > 0.6 GeV/c, to avoid large backgrounds 
at lower momenta. Signal electrons are mostly from pri- 
mary B decays if they are accompanied by a tag electron 
of opposite charge (unlike-sign) . Those with a tag of the 
same charge (like-sign) originate predominantly from sec- 
ondary decays of charm particles produced in the decay 
of the other B meson. Inversion of this charge correlation 
due to B"B° mixing is treated explicitly, and unlike-sign 
pairs with both electrons originating from the same B 
meson are isolated kinematically. With a small model- 
dependence on the estimated fraction of primary elec- 
trons below p* — 0.6GeV/c, we infer the semileptonic B 
branching fraction from the background corrected ratio 
of unlike-sign electron pairs to tag electrons. 

This measurement is based on data recorded in the 
year 2000 with the BABAR detector J| at the PEP-II en- 
ergy asymmetric e + e~ storage ring || at SLAC. The de- 
tector consists of a five- layer silicon vertex tracker (SVT), 
a 40-layer drift chamber (DCH), a detector of internally- 
reflected Cherenkov light (DIRC), and an electromag- 
netic calorimeter (EMC) all embedded in a solenoidal 
magnetic field of 1.5 T and surrounded by an instru- 
mented flux return (IFR). To ensure the high quality of 
the data, we have selected the largest contiguous block 
of events with identical and stable detector conditions 
in the year 2000, corresponding to an integrated lumi- 
nosity of 4.1ffcT 1 collected at the T(45) resonance, and 
0.97 fb" 1 recorded about 40MeV below the T(45) peak 
(off- resonance) . 

Multihadron events are selected by requiring a charged 
track multiplicity of iV c h > 4, or A c h = 4 plus at least 2 
neutral energy deposits above 80 MeV in the EMC. Track 
pairs from converted photons are not included in iV c h, but 
count as one neutral particle. For further suppression of 
non-SB events, we require i?2 < 0.6, where i?2 is the 
ratio of Fox- Wolfram moments H2/H0 |Q. 

The electron momentum measurement and identifica- 
tion are critical for this analysis. For electron candidates 
we require hits in at least 12 DCH layers, and a polar 



angle 9 within the EMC acceptance, i.e. —0.72 < cos0 < 
0.92. To reduce the contamination from photon conver- 
sions and beam-gas background we require the track im- 
pact parameters in the plane perpendicular to the beams 
and along the detector axis to be less than 0.25 cm and 
3.0 cm, respectively. 

The track finding efficiency e t rk is determined from 
data as a function of charged multiplicity, transverse mo- 
mentum, polar and azimuthal angle. For signal elec- 
trons with p* > 0.6GeV/c, the average efficiency is 
(97.1 ±1.1)%. 

Electron identification is based on the ratio of the en- 
ergy in the EMC and the track momentum, -Eemc /p, the 
shower shape in the EMC, the specific energy loss dE/dx 
in the DCH, and the number of Cherenkov photons and 
the Cherenkov angle measured in the DIRC. Muons are 
eliminated on the basis of dE/dx and -Eemc/p- Taking 
into account the correlations between deposited energy 
and shape in hadronic showers, we combine probability 
density functions derived from data samples for each dis- 
criminating variable to construct the likelihood function 
£(£), £ € {e, 7r, K, p}. A track is identified as an electron 
if 



L(e) + 5 L(tt) + L{K) + 0.1 L(p) ' ' 

The weights roughly reflect the relative abundances, their 
exact values not being crucial for electron identification. 

We measure the electron identification efficiency as a 
function of p* and center-of-mass polar angle 9* using 
radiative Bhabha events. For momenta p* > 0.6GeV/c, 
the average efficiency is 92% (see Figure la). However, 
Monte Carlo simulations indicate that relative to radia- 
tive Bhabha events, the identification efficiency in BB 
events is reduced between (4±2)% at low momenta (p* < 
1 GeV/c) and (2±1)% above p* = 1.6 GeV/c. We correct 
the measured efficiency for this momentum-dependent 
difference. 

The misidentification rates for pions, kaons, and pro- 
tons are extracted from control samples selected from 
data. Figure [l]b shows the misidentification probabilities 
rjh per hadronic track, where the relative abundance of 
pions, kaons, and protons is taken from BB Monte Carlo 
simulation. The DCH and DIRC contribute significantly 
at low momenta, while the performance of the EMC in- 
creases with p* . This leads to a minimum of 0.05% for 
r\h at 1 < p* < 1.3GeV/c. The relative systematic error 
is estimated to be 15% from the purities of the control 
samples and the uncertainties in the relative abundances. 

The branching fraction analysis makes use of three 
samples: (1) the tag electrons, (2) unlike-sign and (3) 
like-sign pairs of a tag and a signal electron candi- 
date. Misidentified hadrons and electrons from non- 
Si? (continuum) events, photon conversions, Tr a ,r) — > 
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FIG. 1: Electron identification efficiency e e as obtained from 
radiative Bhabha events (a) and hadron misidentification rate 
r/h (b) as a function of p*. 



7e + e~( "Dalitz") and J/0 , tp(2S) — * e + e~ decays con- 
tribute to the background in all three samples. The 
unlike-sign sample also contains pairs of primary and sec- 
ondary electrons from the same B meson decay. Further 
contaminations to the like- and unlike-sign samples arise 
from decays of r leptons and charmed mesons produced 
in b — > ccs decays. Apart from the correction for unlike- 
sign electron pairs from the same B, which is performed 
in bins of p* only, all background corrections are per- 
formed in bins of p* and polar angle 9*. 

The continuum background is subtracted from all three 
samples. It is obtained by normalizing the observed off- 
resonance spectra by the ratio of on- to off-resonance 
integrated luminosities. The relative systematic error in 
this ratio is estimated to be 0.5%, attributed to variations 
in the detector performance over time. The continuum 
momenta are scaled by the ratio s on j \fs ff to compen- 
sate for the 0.4% lower center-of-mass energy. 

Electrons from photon conversions and Dalitz de- 
cays are identified by pairing them with any oppositely 
charged track with transverse momentum p t > O.lGeV/c. 
We distinguish the two sources of pairs by the distance 
Rpair of the pair vertex from the detector axis. Photon 
conversions are identified by requiring R pa ir > 1.6 cm, a 
pair invariant mass M ee < 100MeV/c 2 , and the trans- 
verse and longitudinal distances between the two tracks 
at the point of closest approach A xy < 0.3 cm and 
A z < 1.0 cm. For Dalitz pairs, we require R pa ir < 1.6 cm, 
M ee < 200MeV/c 2 , A xy < 0.2 cm and A z < 1.0 cm. The 
momentum- and polar angle-dependent pair finding effi- 
ciency, which is obtained from a full detector simulation, 
is low since, in most cases, the momentum of the sec- 
ond track is too small to produce a track in the DCH. It 



varies between 30% and 40% for photon conversions and 
between 20% and 30% for Dalitz pairs. From a detailed 
comparison between data and simulation, including the 
energy spectra of the pairs, the relative systematic un- 
certainties are estimated to be 13% and 19% for the con- 
version and Dalitz background rates, respectively. 

In the unlike-sign sample, electrons from primary and 
charm decays of the same B tend to be produced in op- 
posite directions. Defining p* as the center of the signal 
electron momentum bin, this background is reduced by 
a factor of 24 by imposing the condition 

cos a > 1.0 -p*/ (GeV/c) and cos a > -0.2 (1) 

on the opening angle a of e + e _ pairs, measured in the 
T(AS) frame. Since B mesons are nearly at rest in this 
frame, there is no angular correlation between two elec- 
trons from different B mesons, and the loss in signal 
efficiency can be calculated on the basis of geometrical 
acceptance. 

This selection also eliminates most e + e~ pairs from 
inclusive B — > J/ipX decays. Electron candidates that 
can be combined with an oppositely charged electron to 
form an invariant mass consistent with the J/ip hypoth- 
esis, 2.90 < M ee < 3.15GeV/c 2 , are excluded from the 
tag sample if cos a < —0.2. 

The contribution of unlike-sign pairs from the same B 
decay satisfying Eq. |l| is approximately 2%. After sub- 
traction of background contributions from continuum, 
photon conversions and Dalitz decays, the observed open- 
ing angle distribution (without the requirement) contains 
a flat contribution from electron pairs from different B 
mesons and a contribution from electron pairs from the 
same B, which peaks at cos a = —1. The shape of the 
non-flat background is taken from Monte Carlo simula- 
tion and the relative normalization of the two contribu- 
tions is determined by a fit to the data, which is per- 
formed separately for each 100 MeV/c-wide momentum 
bin below 1.2 GeV/c. The integral over the fitted non-flat 
contribution between the minimal allowed value of cos a 
and 1 is taken as the residual background (Figure |^). 
The very small background above 1.2 GeV/c (0.8% of 
the total contribution) is determined from Monte Carlo 
simulation with a relative uncertainty of 50%. 

We have studied systematic uncertainties in the pre- 
dicted opening angle distributions by varying the branch- 
ing fractions of B — > Dev, B — > D*ev, B — > D**ev 
and non-resonant B — > D^*>ixev decays by one standard 
deviation around current average values [0. Based on 
detailed studies and variations of the fit, the combined 
systematic error for this background is estimated to be 
5%. 

Figure || shows the observed momentum spectra and 
the individual background contributions discussed so far, 
corrected for tracking efficiency; a summary of yields is 
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FIG. 2: Distribution of the cosine of the opening angle of 
unlike-sign pairs for 0.7 < p* < 0.8GeV/c. The points repre- 
sent the data and the histogram is the result of a fit. The 
shaded area represents the estimated contribution of back- 
ground electrons, and the vertical dashed line indicates the 
requirement on the opening angle. 



given in Table |. Following this initial set of background 
corrections, the electron yield is corrected for electron 
identification efficiency. 
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FIG. 3: Total measured spectrum (points) and estimated 
backgrounds (histograms) for signal electron candidates in (a) 
the e + e~ sample, and (b) the e e sample. 



Background contributions from B 



D 



ev e Y decays and B 



DD {S) X, 



e decays are esti- 



mated by Monte Carlo simulation, using the currently 
known branching fractions. Combining B(D S — > Xev) — 
(8.12 ± 0.68)%, which is computed from the average D 



branching fraction B(D a ' + — > Xev) [0 and the lifetime 
ratios t D o,+ /t Ds , with B(B -> D S X) = (9.8 ± 3.7)% §] 
yields B(B -> D s -> e) = (0.80 ± 0.31)%. We take 
the inclusive branching fraction B(B — » DD^X) to be 
(8.2 ± 1.3)% [||. Assuming equal production rates of D 
and D* , but allowing for any ratio in the systematic er- 
ror, we arrive at B{B -> D -► e) = (0.84 ± 0.21)%. To 
estimate the contribution of electrons from r decays, we 
use B(B -> Xtv) = (2.6 ± 0.2)%, B(D S tv) = (5.79 ± 
2.00)% @ and B(t -> e v e u T ) = (17.83 ± 0.06)% @. 
This leads to B(B -> r -> e) = (0.565 ± 0.063)%, where 
the t lepton originates either directly from a B decay 
or a B — > D s — > r cascade. The background from J/t/i 
and tp(2S) decays into two electrons is also estimated 
from Monte Carlo simulation, with B(B — > J/'i/i — > e + e~) 
= (6.82 ± 0.38) x 10~ 4 and 6(5 -> ^>(2S) -> e+e") = 
(3.1 ±0.6) x 10~ 5 §. 

The tag electron sample is first corrected for contin- 
uum background and hadron misidentification. The re- 
maining background is from secondary decays of charm 
particles and unvetoed Jji\) — + e + e~ decays. All these 
contributions are estimated by Monte Carlo simulation, 
leading to the background-subtracted number of tag elec- 
trons N Ug = 304,048 ± 880(stat) ± 2,100(syst) (Table |), 
including a correction for signal loss due to the J/tp -veto. 

Due to B°B° flavor oscillations, electrons from primary 
B decays and B — > DX, D — * e~v e Y cascades contribute 
to both unlike- and like-sign spectra. Denoting the effi- 
ciency of the opening angle cut as e a (p*), their p* distri- 
butions can be written as 



1 dN+- 
e a (p*) dp* 



dX B ^Xev 

dp* 

dN B _Xev 



(i-x) + 



dN 



B^D-^Xeu 



dp* 



dp* 



dp* 



where x is the product of the B°B° mixing parameter 
Xo = 0.174±0.009 § and f = B(Y(AS) -► B°B°). Since 
the measured ratio of charged to neutral T(4S) decays is 
consistent with unity [jio), we assume /q = 0.500 ±0.025, 
where the error is taken from |ic[| . We use these lin- 
ear equations to determine the primary electron spec- 
trum from B decays, dNs^Xev I dp* . Integration of this 
spectrum between 0.6 and 2.5 GeV/c yields Ns^xev = 
25,070 ± 410(stat). Using Monte Carlo simulation, wc 
determine the relative efficiency for selecting events with 
two electrons compared to events with a single tag to be 
e e vt = (98.0 ± 0.5)%. Together with the polar angle ac- 
ceptance e geom — 84%, we obtain the partial branching 
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TABLE I: Electron yield for the three samples and corrections with statistical and systematic errors. 





(1) tag sample 


(2) e + e sample , cut on a 


(3) e ± e ± sample, all a 




1.4 < p* < 2.3 GeV/c 


0.6 < p* < 2.5 GeV/c 


0.6 < p* < 2.5 GeV/c 


On T(4S) 


395,791 ± 630 


14,692 ±120 


10,838 ±110 


Continuum 


82,073 ± 590 ± 410 


1,301 ±76 ± 7 


939 ± 64 ± 5 


_i_ 

7 — > e T e 


561 ±23 ± 140 


283 ± 40 ± 37 


856 ±82 ±110 


o -t- — 
r/, 7r — ► 7e e 


92 ± 9 ± 23 


51 ±22 ±10 


80 ±82 ±15 


Faked e 


1,455 ± 140 ± 360 


136 ± 16 ± 20 


348 ± 48 ± 52 


e from same B 




317±7±16 




Yield before and 


311,610 ±870 ±570 


12,603 ±150 ±46 


8,616 ±180 ±120 


after eff. corr. 




14,134 ± 180 ± 170 


9,734 ± 190 ±200 


B — > r — > e 




353 ±17 ±42 


93±9±11 


r> — > _L7 S — > e 




zn Lu zn 11U 


79 -4- Q -4- 98 


B -> D ->e 




226 ± 16 ± 57 


65±8±16 


Secondary tags 


8,073 ±91 ±2, 000 


296 ± 17 ± 74 


886 ± 29 ± 220 


e from J/i/> or ip(2S) 


1,925 ±42 ±120 


77±8±5 


119±10±7 


e removed by J/ip veto 


-(2,435 ± 50 ± 220) 






Net e yield 


304,048 ±880 ±2,100 


12,890 ±180 ±230 


8,500 ±200 ±300 



fraction 



B[B -> Xev, p* > 0.6 GeV/c) 



Xev 



Ntag £brem ^evt £geom 

= (10.24 ± 0.17(stat) ± 0.26(syst))% , 



which includes a correction for the small loss of electrons 
due to bremsstrahlung in the detector material and the 
limited momentum resolution, l~et, rem — (2.20±0.35)%. 
The contributions to the systematic error are listed in 
Table || 

Figure V& shows the momentum spectrum of 
primary electrons. 
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FIG. 4: Momentum spectrum of electrons from decays 
B — » Xev after correction for efficiencies and external 
bremsstrahlung, with combined statistical and systematic er- 
rors. The curve indicates the fit used for the extrapolation to 
p* = 0. 



To determine the total semileptonic branching fraction, 
we need to extrapolate the spectrum to p* = 0. This is 
achieved by fitting the data to the sum of the spectra 
from the various exclusive decays. We use a parameter- 
ization of HQET-derived form factors jH], [l2| to model 
the decays B — > Dev and B — ► D*ev, and the work of 
Goity and Roberts jl3| for non-resonant B — * D^nev 
decays. Semileptonic B decays to D**ev and charmless 
mesons are described by the ISGW2 model jjj, which is 
also used as an alternative description for the processes 
B — > Dev and B — > D*ev. Photon radiation in the final 
state is modeled by PHOTOS [|l5). The relative contri- 
butions of the different exclusive decay modes are con- 
strained to be within two standard deviations of the mea- 
sured average branching fractions jjj . The best estimate 
for the extrapolation factor is 1 ± n = 1.061 ± 0.009, 
where the error accounts for the observed variations of 
the fit results for different decay models and branching 
fractions. This extrapolation leads to a total semileptonic 
branching fraction Bsl of 

B(B -► Xev) = (10.87 ± 0.18(stat) ± 0.30(syst))% . 



One of the limiting factors of this analysis is the back- 
ground at low momenta, especially semileptonic decays of 
charmed mes ons produced in b — > ties decays. As shown 
raising the minimum momentum require- 



in Table |TII 

ment p* min reduces the systematic uncertainty due to this 
background substantially, but also increases the error on 
the extrapolation to p* = 0. We choose p* min — 0.6 GeV/c 
for the final result, since the systematic error is com- 
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TABLE II: Impact of systematic uncertainties on Bsl- 



Source 


AB S l(%) 


Source 


AB(%) 


e efficiency 


0.144 


B — > t — > e 


0.044 


B -> D a -> e 


0.130 


7 — » e + e~ 


0.042 




0.120 


^■bT G7TL 


0.039 


extrapolation 


0.092 


faked e 


0.024 




0.075 


e from same B 


0.022 


fl -» D -» e 


0.067 


7T°, 77 — > 7 e + e~ 


0.014 


mistagged e 


0.061 


continuum 


0.008 


/0X0 


0.059 


J/0,i/>(2S) -> e+e" 


0.003 




0.054 






Total 






0.296 



parable with higher values of p* 
dependence is significantly lower. 



while the modcl- 



TABLE III: Determination of k, Bsl, and the contributions 
to the systematic error for different signal electron momentum 
cut-offs. All numbers are stated in percent. 



?w[GeV/c] 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


n 


3.8 


6.1 


9.3 


13.6 


19.2 


27.2 


Bsl 


10.79 


10.87 


10.87 


10.82 


10.80 


10.93 


ASsl(7,7t°) 


0.07 


0.04 


0.03 


0.02 


0.02 


0.01 


ABsL^trk) 


0.12 


0.12 


0.12 


0.12 


0.12 


0.12 


AB SL (e eff.) 


0.15 


0.14 


0.14 


0.12 


0.11 


0.10 


ABs L (B->D a ) 


0.17 


0.13 


0.09 


0.06 


0.05 


0.04 


AB S l{B->D) 


0.10 


0.07 


0.05 


0.03 


0.02 


0.01 


ABsl(B^t) 


0.05 


0.04 


0.04 


0.03 


0.03 


0.02 


ABs l (extrapolation) 


0.06 


0.09 


0.13 


0.19 


0.25 


0.33 


ABsl (other) 


0.15 


0.14 


0.14 


0.15 


0.15 


0.17 


ASsL(syst) 


0.33 


0.30 


0.29 


0.30 


0.34 


0.41 


AZ3 S L(stat) 


0.21 


0.18 


0.16 


0.16 


0.15 


0.15 



Based on the work by Hoang et al. 16 , we relate the 
decay rate and the modulus of the CKM matrix element 
V cb by 

\V cb \ =(41.9 ±2.0) x 1(T 3 



x y/B(B X c ev)/ 0.105 y/l.6ps/T B . 

Using t b = (1.601±0.021) ps and B(B -> X u ev) = (1.7± 
0.6) x 10- 3 § , we obtain \V cb \ = 0.0423 ±0.0007(exp)± 
0.0020(theory). 

In conclusion, we have used electrons in T(4S) de- 
cays tagged by a high momentum electron to measure 
B(B -► Xev)= (10.87 ± 0.18(stat) ± 0.30(syst))%. This 
measurement is largely model- independent. The result 
is in agreement with previous measurements m, O , but 



the systematic uncertainties are reduced. However, the 
poorly known branching fractions in B and Dr s \ decays 
lead to significant systematic uncertainties in the back- 
ground subtraction. The resulting measurement of \V c b\ 
remains dominated by theoretical uncertainties. It has 
recently been shown that non-perturbative effects can be 
assessed by measurements of moments of inclusive distri- 
butions pl| . 
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